Modelling vegetation land fragmentation in urban areas of Western Province, Sri Lanka using an Artificial Intelligence-based simulation technique

Vegetation land fragmentation has had numerous negative repercussions on sustainable development around the world. Urban planners are currently avidly investigating vegetation land fragmentation due to its effects on sustainable development. The literature has identified a research gap in the development of Artificial Intelligence [AI]-based models to simulate vegetation land fragmentation in urban contexts with multiple affecting elements. As a result, the primary aim of this research is to create an AI-based simulation framework to simulate vegetation land fragmentation in metropolitan settings. The main objective is to use non-linear analysis to identify the factors that contribute to vegetation land fragmentation. The proposed methodology is applied for Western Province, Sri Lanka. Accessibility growth, initial vegetation large patch size, initial vegetation land fragmentation, initial built-up land fragmentation, initial vegetation shape irregularity, initial vegetation circularity, initial building density, and initial vegetation patch association are the main variables used to frame the model among the 20 variables related to patches, corridors, matrix and other. This study created a feed-forward Artificial Neural Network [ANN] using R statistical software to analyze non-linear interactions and their magnitudes. The study likewise utilized WEKA software to create a Decision Tree [DT] modeling framework to explain the effect of variables. According to the ANN olden algorithm, accessibility growth has the maximum importance level [44] between -50 and 50, while DT reveals accessibility growth as the root of the Level of Vegetation Land Fragmentation [LVLF]. Small, irregular, and dispersed vegetation patches are especially vulnerable to fragmentation. As a result, study contributes detech and managing vegetation land fragmentation patterns in urban environments, while opening up vegetation land fragmentation research topics to AI applications.


Introduction and literature review
According to the findings of Abdullah et. al. [41] the five most populous cities in Bangladesh-Dhaka, Chattogram, Khulna, Rajshahi, and Sylhet-have lost more than USD 628 million as a result of specific ecosystem degradation, with the decline of green and blue land areas being the primary cause of this loss. The urban heat island effect is mostly caused by the degradation of vegetated land. For example [42], study demonstrates that due to human limits on agricultural operations, enforcing lockdown following the COVID-19 has increased the amount of vegetative land cover in the Indo-Gangetic basin, India, throughout the winter. This increases the amount of evapotranspiration, which accounts for about half of the region's cooling [42]. Additionally, [43] shows that in five main cities in Bangladesh, including Dhaka, Chittagong, Khulna, Rajshahi, and Sylhet, a lack of green space and a high proportion of impermeable surface are significant contributors to urban heat. Furthermore, vegetation works as a natural flood barrier; hence, when it vanishes or changes as a result of development, there are serious flooding problems in urban areas. According to [44] research, mangroves, marshes, and coastal vegetation are essential for preventing coastal flooding in the world. The irregularity and complexity of fragmented lands can affect land-use efficiency, lowering the economic return on agricultural land [3,45,46]. Also, the urban structure's fragmentation might operate as a barrier to social connection [47]. As a result, vegetation land fragmentation has an impact on sustainability because of the physical structure of land use [1,11]. Therefore, urban planners are now more interested in exploring vegetation land fragmentation alongside landscape and environmental planners [5]. So far, national, regional, and local level empirical studies have been used to study vegetation land fragmentation, while noting the extremely limited yet successful research efforts on vegetation land fragmentation [48]. However, there is an unmet demand for multi-factor modeling of vegetation land fragmentation in urban environments [2,33]. Because vegetation land fragmentation is caused by various variables in a complex environment, the traditional linear regression model is unable to convey the relevance of influencing factors [1,21]. To examine nonlinear interactions in the urban setting, several recent research has used AI-based modeling tools [1,22,[47][48][49]. However, developing an AIbased modeling framework with multiple factors to explicitly model vegetation land fragmentation is a knowledge gap that needs to be investigated to control vegetation land fragmentation and its consequences [3,27]. Even though numerous researchers have discovered factors that contribute to vegetation land fragmentation, it is critical to understand the types of relationships that exist between vegetation land fragmentation and influencing factors (rules) to make spatial planning judgments. Most of the recent research have identified the variables, but there has been little focus on quantitatively explaining the relationship [50]. To develop efficient land-use regulations and guidelines, it is necessary to determine the magnitude of variables' impacts on vegetation land fragmentation [27,47]. even though successful research attempts on vegetation land fragmentation have focused on forest land cover rather than the entire vegetation land cover since they are more oriented to the ecological domain [51]. In addition, vegetation land fragmentation occurs in both rural and urban areas while generating the same kind of impacts [32,36]. Since the study focuses on the spatial pattern of vegetation land use rather than the ecological process, this study focuses more on vegetation land fragmentation in urban areas.
Sri Lanka is a South Asian country with a diverse range of green and blue land uses; and has become a gateway to Asia as a result of regional development projects-such as Port City Colombo, Hambantota Harbor, and Airport Development. Rapid infrastructure development projects, such as expressway construction and regional infrastructure development, have resulted in a reduction in the percentage of land covered by vegetation in Sri Lanka [52]. Western Province in Sri Lanka has seen significant urbanization and has lost most of its natural cover [14,53]. Various research has used GIS and remote sensing techniques such as NDVI [54][55][56], Satellite image classifications, and spatial metrics to explore land cover change in Sri Lanka [57,58]. Some researchers have even used AI-based applications such as the SLEUTH model at the regional level to anticipate land cover changes [59]. However, a handful of research have looked into vegetation fragmentation in urban areas, in Sri Lanka. In addition, one study has focused solely on paddy field fragmentation and its economic consequences [46]. Because Sri Lanka is rich in biodiversity and forest, wetlands and other vegetation covers, it is important to investigate the level of vegetation fragmentation. Due to the growing urbanization patterns in the Western Region, this research focuses on vegetation fragmentation in the Western Province; because studies of vegetation fragmentation are necessary for determining implications for reducing fragmentation's negative effects and ensuring sustainable development in urban areas [4,27,60]. Urban planning agencies such as the Urban Development Authority, the National Physical Planning Department, and the Central Environmental Authority, as well as local governments, require tools and model applications to forecast future scenarios in vegetation fragmentation and identify their factors in order to develop effective policies to control the level of fragmentation in Sri Lanka and to make society more sustainable. Although there are already developed AI-based applications such as SLEUTH, MOLUSCE, they are oriented to analyze the urban growth or land cover changes. Therefore, it is essential to develop an AI-based model framework to model vegetation land fragmentation not only in Sri Lanka but also for other countries.
The objective of this study is to develop an AI-based simulation framework to simulate vegetation land fragmentation in urban environments. The non-linear relationships and magnitude of impacting elements were identified using a feed-forward ANN in this study. Under the AI modeling technique, the study also built a DT modeling framework to explain the effect of many variables on vegetation land fragmentation. The study exclusively used data from Sri Lanka's Western Province to create the model. As a result, the study's scope is restricted to the Western Province. This work, on the other hand, might be viewed as the first attempt to measure and simulate vegetation fragmentation in the Sri Lankan environment. Furthermore, the methodological scientific relevance and contribution of AI-based modeling approaches adds to the incentive to this paper.

Study area and data sources
Western Province, Sri Lanka, was chosen as the case study considering data availability and the capacity to meet the sample size requirements of AI models. According to the 2012 census and statistics report in Sri Lanka, the Western Province covers 3684 km 2 and has a population of 5.85 million people. It has the densest urban population and covers 46 percent of the total built-up area. With significant urban and infrastructural developmental projects such as expressways, port and airport extensions, and power stations, the Western Region serves as Sri Lanka's commercial hub. It also contains 2505 Grama Niladhari Divisions (GND), which are local administrative borders, and these 2505 GNDs were utilized as the AI model's sample in this study.
For the study, primary and secondary data sources were employed. The vegetation and built-up land-uses in the Western Province were mostly extracted from USGS Lands at 7 satellite images in 2000 and 2010. The elements of vegetation land fragmentation were calculated using data from the Survey Department's Road network and contour layers, as well as data from the Census and Statistic Department of Sri Lanka's population and building units. Table 1 shows the data sources for each study variable in further detail.

Method
The study is based on a review of the literature, which included both theoretical and empirical investigations, to determine the factors that cause vegetation land fragmentation from the beginning [Fig 2]. The identified factors were then conceptualized as patch, corridor, matrix, and other factors based on the theoretical explanation of landscape ecology. Fig 3 illustrates the conceptual framework for vegetation land fragmentation in this study. The study relied on 20 different parameters, which are included in Table 1, and used USGS satellite images of Western Province in 2000 and 2010 to classify them into vegetation and built-up classifications (with over 75% classification accuracy). At the GND level, each factor is calculated using geoprocessing tools available in ArcMap and QGIS. Even though the components were conceived as patch, corridor, matrix, and other categories in the study, the geographic entity employed was GND, Sri Lanka's smallest administrative entity. The availability of secondary data was also taken into account in this investigation. As a result, the model outputs are also tied to GND administrative borders rather than spatial features such as patches, corridors, or matrix. The vegetation fragmentation was quantified using FRAGSTATS (Landscape Division Index). During the database development, existing local reports determined some of the quantitative parameters. Before defining the models, Principal Components Analysis was used to exclude the multi-correlated and least-correlated factors. After all, the study used 1750 training data and 750 testing data to design and validate the ANN and DT models. Finally, using the ANN, it calculated the future LVLF in Western Province. The next sections will provide extensive explanations for each process.
The landscape division index was used to calculate the LVLF at the class level. In the FRAG-STATS software, it is one of the landscape metrics. By detecting random pixels in each area, the landscape division index assesses the division of the same land use category. It determines the likelihood of two random pixels not being in the same patch [64].
The complete landscape area is denoted by the letter A [61]. aij is the area of patch ij [61].   The landscape division index has a value ranging from 0 to 1 [61]. The value 0 indicates that the area is not fragmented. When the value is close to 1, the land is highly fragmented. The value will be 0 if j and i pixels are in the same patch [61]. If the pixels j and i are in distinct patches, the value ranges from 0 to 1 [61]. The division value will be closer to 1 if i and j pixels are located inside smaller areas.

AI-based modeling frameworks.
The complex urban environment processes and challenges have been modeled using AI-based modeling frameworks [65]. The bulk of research attempts to do away with the Linear Regression model, since it is insufficient when dealing with a complex environmental issue [66]. The advantages of using AI-based modeling frameworks include the ability to model with several variables and the ability to identify the cluster influence of the variables [67]. Consequently, AI models can be used to identify nonlinear connections. To meet the research objectives, this study used both ANN and DT.

Artificial Neural Network [ANN].
In this study, the first AI-based approach was ANN, which was used to predict vegetation fragmentation. The approach of an intelligent system that imitates the behavior of the human brain [neurons] is known as ANN [68]. The first mathematical neural model was established in 1943 by Warren McCulloch and Walter Pitts [69]. Input, hidden, and output layers [70]are the three primary layers of an ANN, as shown in Fig 5. Fig 6 represents a rudimentary neural network and its output-generating function. The input variables are X1, X2, and X3, and the weights of respective inputs are W1, W2, and W3. The bias node is 'b,' and the output is 'Y.'. The activation function, which is the neural network's mapping process, is utilized to smooth the output results [67]. As a result, the output function is the multiplication of weight, input, and bias [69].
There are two types of neural networks: supervised and unsupervised [70]. This research uses a supervised neural network to train the algorithm from provided outputs and inputs by changing the predicted output to the actual output while updating the weights according to neural network rules. Data categorization, change detection, grouping, forecasting, and predictions are all possible using supervised neural networks (Deep Learning) [71]. It uses hidden layers to train the existing nonlinearity between input and output variables and predicts the output [72]. Therefore, it can determine the level of significance of input variables in forecasts [68]. It is more useful for simulating spatial dynamics and identifying changes in variables [73].
Before the training process, the data was standardized using the lowest and maximum values. x is the original data value, and y is the matching data value after the maximum and lowest values of x have been normalized [74]. The major packages used in R Studio for ANN preparation are Keras, Neuralnet, and Tensorflow. The default algorithm of the R language's rprop + algorithm, which is based on resilient backpropagation, was utilized in the study. To identify the minimal error function, it changed the weights of input nodes [70]. Because the neural network's linear output function is false, the default logistic activation function was used, with output values ranging from 0 to 1. Mean Square Error [MSE] [75] and Root Mean Square Error [RMSE] [76], which examine the error between predicted and tested data, were used to evaluate the model [70].
Using the R neuralnet package's Olden and Lekprofile algorithms, the study determined the feature significance of variables and rules. Olden is an algorithm that uses the raw value of input-hidden and connection weights of hidden outputs of each input and neuron to indicate the importance of variables [77]. The total of all hidden neurons is then calculated. Each variable's negative and positive importance levels can be displayed by the olden algorithm [78]. In the developed ANN, the lekprofile method is used to verify the sensitivity of each variable to the independent variable [68]. Other factors are kept constant while determining a variable's sensitivity. For each variable, the Lekprofile algorithm generates a graph with probable sensitivities or relationship curves. The complete technique of framing the ANN model through R is shown in Fig 7. 2

.2.2.2. Decision Tree (DT).
In machine learning, the DT is a strong classification and modeling technique [51] which is employed as the second AI-based technique in this study. It builds a tree structure similar to a flow chart and is classified as supervised machine learning [79]. Individuals' behavior or decisions about one another can be predicted by the DT [80]. The DT was calculated using WEKA software in this investigation. The technique utilized was the J48 algorithm, which solely classified the model's significant factors. The data were first normalized using the equation of minimum and maximum, and then the arff. file was generated.  [82] were used to validate the DT model. Fig 8 describes the formulation and validation procedure for DT models.
P agree is about the observed agreement [81]. P chance is agreement expected by chance alone [81]. More than 70% kappa statistic value would represent the validated model in terms of higher accuracy [81].

RAE ¼ j
Observed Actual Value À Expected value Expected Value j:100% ð8Þ RAE examines the absolute error to actual value ratio [82]. The discrepancy between the actual value and the measured value is known as absolute error [82]. The model is considered accurate if the RAE is less than 20%.

RRSE ¼
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The difference between the predicted (Pij) and the target value is also calculated by RRSE (Tj) [82]. It divides the total squared error of the simple predictor to normalize the total squared error.

Results
The study first measured the LVLF and factors of LVLF, then framed the ANN and DT models to discover the variables of LVLF and their relationships, as indicated in the method. Finally, the study used ANN to simulate future LVLF in Sri Lanka's Western Province from 2010 to 2030. Therefore, the results of the ANN and DT models, as well as the simulation results, will be explained in the following section.

ANN
The architecture of the neural network is depicted in Fig 9. Accessibility Growth (AG), Initial Vegetation Large Patch Size (IVLPS), Initial Vegetation Land Fragmentation (IVLF), Initial

PLOS ONE
AG is the positively most significant factor of the LVLF, according to the olden algorithm of ANN (Fig 10). Between the ranges of 50 and -50, the significance level of AG is 45. The IVPA (distance between patches) is the second-most important positive factor which is 23. IVLPS has the greatest negative importance level -14, whereas IBD has the second-highest negative importance level which is -11. The remaining components have negligible positive significance ratings. According to the findings, AG is the most important factor in corridor considerations, and it is also the highest. IVPA and IVLPS are two patch-related parameters that are significantly important.

PLOS ONE
The rules of each explanatory variable may be determined using ANN's 'lekprofile' analysis. Fig 11 illustrates a composite representation of the LVLF's rules. AG interferes favorably with LVLF, whereas increasing AG causes LVLF to rise dramatically. The relationship between IBD and LVLF is negative, indicating that the two variables do not move in the same direction. With increasing building density, LVLF eventually diminishes. The IBLF and the LVLF have a positive relationship. The LVLF gradually increases as the LBLF grows. When the shape of the vegetation becomes linear patches, LVLF continuously rises, indicating that IVC and LVLF have a positive connection. The IVSI and the LVLF have a positive correlation. LVLF steadily rises as the shape of the vegetation patches becomes more uneven. If the relationship between IVLF and LVLF is positive and LVLF is rising in lockstep with IVLF, it is a positive indicator. IVLPS has a negative relationship with LVLF, and LVLF gradually diminishes as the size of the

PLOS ONE
vegetation patch grows larger. The LVLF and the IVPA (distance between patches) have a positive relationship, meaning that increasing the distance between vegetation patches gradually raises the LVLF.

Decision Tree (DT)
The results of the DT model can be used to determine the scenarios of LVLF (Fig 12) and its factors. The DT analyses the link between variables and various scenarios of LVLF with the significant factors as discovered factors in the neural network. The ANN model's greatest priority level factor is AG and the DT model double-proves AG as the root of LVLF. It explains why, in each location, LVLF is particularly high when AG is very high. If the AG is low and the IBD is high, the LVLF is mild. When AG is low and IBD is high, the model uses the IVPA (distance between patches) variable as a branch to explain the various scenarios. If IVPA is high in comparison to earlier connections, LVLF will be moderate. With IVLPS linking to the extremely low AG, high IBD, and very high IVPA, the model has formed a new branch. In keeping with the foregoing relationships, if IVLPS is very high, the LVLF is very low. With the IVLF variable and the IVC variable, another alternative possibility has emerged. In keeping with the connections, if the IVLF is very high, the LVLF is also quite high. With the link of preceding factors, if IVC is very high, LVLF is also very high. Therefore, the DT is a model that can be used to

PLOS ONE
determine LVLF using these variables' connections (factors). The numeric values of each variable related to very low to very high are shown in Table 2.
The results of the model validation demonstrate that 98 percent of the data were properly categorized. RAE and RRSE are 12 and 29 percent respectively. As a result, the model has a greater level of accuracy in terms of validation methods. Fig 13 indicates the percentage of accuracy levels in the tested model using training data. It demonstrates that a minimum of 1000 observations is required to anticipate accurate results. The model accuracy improves as the number of training observations increases from 1000 to 1750. The confusion matrix (Table 3) shows how effectively the DT model classifies the various levels of categories inside the model. Significantly, all the cases with very high LVLF have been accurately identified, whilst the remainder of the groups contain slight misclassifications.

Future simulation
In comparison to the other two districts in the Western Province, the future simulation from 2010 to 2030 (Fig 14) shows a large increase in LVLF in the Gampaha District. It also shows the increase in LVLF along Western Province's expressways and their interchanges. In the Colombo core region, LVLF exhibits a downward trend from 2000 to 2010. The Gampaha District and the Colombo outskirts, on the other hand, have seen a major increase in LVLF.

Discussion
The ANN and DT models were used to simulate vegetation land fragmentation and for identifying important factors and non-linear relationships. A supervised feedforward neural network [Deep Learning] was employed to investigate the variables of vegetation land fragmentation and their behavior, model future vegetation land fragmentation. The DT [supervised classification] model was employed in identifying possible scenarios in vegetation land fragmentation.
The findings of the ANN's olden and lekprofile algorithms, as well as the DT's J48 algorithm, were stated to meet the research objectives. According to the ANN olden algorithm, AG verifies the maximum importance level 45 between the ranges of -50 and 50. When compared to the recent research covered in this study, this is a surprising discovery because many studies regard accessibility as a road network [83,84].
Further the study calculated accessibility using the centrality measures of betweenness centrality and closeness centrality, which is one of the most recent ways of assessing accessibility using space syntax [85]. The rapid road development in the suburbs, such as expressway extensions, is the explanation for this discovery. The expansion of the road network improves accessibility in both urban and suburban areas [85,86]. Consequently, the model clearly addresses the increase in accessibility expansion and its beneficial effects on vegetation land fragmentation. Urban sprawl has a direct influence on the division of vegetation areas, the study's findings point to increased accessibility as the primary cause of vegetation fragmentation [36,87]. The findings of this investigation back up the notion of landscape ecology and the fragmentation process as a contribution of this study. With the LVLF, the second-highest positive significance level 23 is associated with vegetation patch association or distance between vegetation patches, which is consistent with the findings of previous research [7,34]. Although earlier research has looked at the configuration factors to quantify the LVLF [3,88,89], this study considers them to be factors of vegetation land fragmentation, as the configuration of the vegetation patches always causes fragmentation. Therefore, the study contributes our understanding of how to account for patch-related factors when modeling the LVLF, which is another contribution to the existing knowledge Patch association, for example, is regarded as one of the drivers of vegetation fragmentation [89], although it is also considered one of the features of vegetation fragmentation in many studies. Further the findings reveal that the association between patches is a significant factor and that reducing the association would raise the LVLF. Patch association is also important in terms of species movements, resource distribution, and ecosystem service, according to basic landscape ecological theory [89,90]. The hypothesis goes on to say that the loss of connection or association will have a detrimental influence on biodiversity, and it encourages nearby land-use types to be converted [29]. As a result, if vegetation patches are closer together, the likelihood of switching vegetation land use is lower than if patches are separated by greater distances [29,90]. However, other factors, such as the size of the vegetation patch, might influence the distance between patches [28,29]. As a result, more research is needed to identify the relationship between patch association with fragmentation and the composition features of the patch. On the other hand, the highest level of negatively affecting factors on the LVLF is the size of the vegetation large patch. The findings are consistent with previous research [40]. Other factors that are positively influencing LVLF based on their non-linearity, such as vegetation shape irregularity, vegetation shape circularity (circle or linear), built-up land fragmentation, and vegetation land fragmentation of the initial year, are similar to fundamental landscape ecology explanations [29,30,32]. The building density of the initial year and LVLF, on the other hand, are not in the same direction as prior research explored, which is consistent with the previous study [50]. The influence of urban sprawl in suburban regions might be the cause of this observation [36,91,92]. In comparison to Colombo core regions, suburban areas in Western Province do not have a higher level of building density [93]. However, because of expressway construction and road network expansions in the Western Province, suburban regions have become more accessible. Therefore, new buildings with lower density are emerging [93]. Owing to new developments, existing vegetation patches in suburbs have been greatly fragmented compared to highly densified areas, according to a previous study [93]. Consequently, the study finds that building density and LVLF have a negative relationship. However, carrying out a local level research for future studies to validate these results can yield a more accurate outcome. Significantly, the research shows that corridor-related parameters are more important than patch and matrix, which opens-up further research opportunities in the field of land fragmentation. Both patch and matrix-related parameters, in contrast, have a considerable impact on LVLF. Specifically, rather than measuring several indices to assess fragmentation, this study identifies vegetation fragmentation and LVLF aids in determining the factors and their non-linear interactions. The analysis solely grouped the components under patch-corridor-matrix with numerical values into GND borders in this case. However, the study does not give exact variables in the form of geographical corridors, patches, or matrixes. Due to this constraint, the model findings for each category cannot be spatially interpreted in this study. Therefore, future research should look into employing exact geographical entities rather than confining itself to the border, as this study did; as it can help to develop a spatial model like Futures, SLEUTH, unlike statistical models. In model applications, the study claims that validation of the ANN displays the strongest level of prediction accuracy, since MSE and RMSE are 0.025 percent and 1.574 percent, respectively. After all, the DT model also shows 12% of RAE and 29% of RRSE. The findings further confirm that similar to environmental studies, the AI model can be used to simulate vegetation land fragmentation [1,2,47]. Furthermore, the future simulation depicts an increasing trend of LVLF along expressways in Sri Lanka's Western Province, demonstrating the relevance of corridor-based variables in promoting LVLF.
The Sri Lankan government has recently placed a greater emphasis on expressway construction, intending to increase regional transportation through huge infrastructure projects [94,95]. In the Western Province, for example, the Southern Expressway, the outer-circular highway, and the Katunayake Expressway have all been completed. The Kandy Expressway is also under development and has begun in the Western Province. Therefore, Western Province serves as the country's transportation center. Consequently, rapid land clearances and divisions associated with expressway constructions may be the primary cause of vegetation land fragmentation in Sri Lanka's Western Province and these findings consist with the study of the effect of social and environmental factors on expressway construction in Sri Lanka [95]. However, the research is confined to the Western Province of Sri Lanka, and it solely looks at vegetation land fragmentation. Thus, it is critical to investigate a variety of case studies and land-use categories such as forest and paddy separately [28,29] while examining fragmentation research, because the fragmentation level may be determined by a variety of factors depending on the land use type [28,29,95]. It is also important to evaluate the case study's various aspects. Furthermore, by referring to local approaches, modeling findings must be validated in the local context. Therefore, future researchers may use the same technique in dissimilar case study areas at the regional, local, and site levels to generalize the study's findings while focusing on land fragmentation in other land use categories other than vegetation. Additionally, future research can investigate the non-linear relationship between vegetation land fragmentation and AI models other than ANN and DT, as well as other factors of vegetation land fragmentation in different categories, in order to identify new determinants of land fragmentation. It is far more necessary to adhere to geographical rather than administrative boundaries because it has the potential to generate several inaccuracies in the analysis findings. Due to the lack of secondary data, this study relied on administrative borders. However, for extremely accurate findings and to construct a spatial model, future researchers need to incorporate spatial data such as patches, corridors, and matrix.
The main contribution of this study is the development of an artificial intelligence-based simulation framework for simulating vegetation land fragmentation in urban regions with acceptable accuracy. The study provides numerical evidence for a nonlinear relationship between factors and land fragmentation caused by vegetation. The study also develops a tree structure to fully describe the phenomena of vegetation land fragmentation. This study contributes methodologically by utilizing AI-based technologies such as ANN and DT to understand complex, non-linear interactions and assess and simulate vegetation land fragmentation.

Conclusion
The study was successful in building an AI-based simulation framework for modeling vegetation land fragmentation in urban environments, as well as in attempting to discover variables that influence vegetation land fragmentation through non-linear connections. Supervised feedforward ANN (Deep Learning) was used to identify the drivers of vegetation land fragmentation, understand how they behave, and model future vegetation land fragmentation. The DT (supervised classification) model was used to identify various possible scenarios in vegetation land fragmentation.
The implications of study's findings are valuable for urban planners who want to learn more about vegetation land fragmentation. Researchers interested in vegetation land fragmentation could utilize this work to consider vegetation land fragmentation from an urban planning viewpoint and model it in a new context. In urban and regional planning, forecasting future land-use changes, identifying urban growth patterns, and assessing the influence on natural vegetation are difficult challenges. As a result, planners may successfully simulate those occurrences using the developed modeling framework; and use the developed DT diagrams and rules to get quantitative insight into vegetation land fragmentation phenomena and regulate them with strategic actions. The study discovered that increasing accessibility has a substantial influence on land fragmentation in vegetation.
Furthermore, fragmentation threatens tiny irregular and scattered vegetation patches. Therefore, the planner can devise measures to reduce vegetation land fragmentation, such as creating interconnected vegetation corridors and reducing accessibility expansion in more sensitive vegetation. Therefore, by examining the context's land uses, planners may propose solutions to prevent the impacts of vegetation land fragmentation. It is difficult to assess the efficiency of spatial planning initiatives in terms of landscape change and fragmentation in contemporary practice. Therefore, the study's findings and method can also be used to evaluate the impacts of existing planning attempts.
Data unavailability and data distortions are major drawbacks of the study. Further, the outcomes of the study cannot be interpreted spatially due to data limitations. Future research can build on this framework in other case study sites with various land use categories. It would also be interesting if future academics could put this system through its paces with diverse AI models. This project will continue the ground verification of the results and the modeling framework. Developing a spatial model, moreover, is vital for successfully interpreting spatial dynamics. Despite limitations, the findings of the study provide a framework for quantifying, analyzing, and modeling vegetation land fragmentation, allowing planners to assess the current situation, forecast future trends, and develop effective spatial planning strategies and land monitoring mechanisms to achieve long-term sustainable development.